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We prove that the quasiclassical tau-function of the multi-support solutions to matrix models, 
proposed recently by Dijkgraaf and Vafa to be related to the Cachazo-Intrilligator-Vafa superpo- 
tentials of the M — 1 supersymmetric Yang-Mills theories, satisfies the Witten-Dijkgraaf-Verlinde- 
Verlinde equations. 



o 

o 
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^ . 1 Introduction 

^ The Witten-Dijkgraaf-Verlinde-Verlinde (WDVV) equations in the most general form can be writ- 

Q ' ten as system of algebraic relations 



TiT-j^Tk = TkT~j^Tu y I,J,K (1) 

for the third derivatives 

of some function ^(T). Have been appeared first in the context of topological string theories PP, 
they were rediscovered later on in much larger class of physical theories where the exact answer for a 
multidimensional theory could be expressed through a single holomorphic function of several complex 
variables H H [71 El . 

Recently, a new example of similar relations between the superpotentials of = 1 supersymmetric 
gauge theories in four dimensions and free energies of matrix models in the planar limit was proposed 
iniEl- It has been realized that superpotentials in some M = 1 four-dimensional Yang- Mills theories 
can be expressed through a single holomorphic function [U] that can be further identified with free 
energy of the multi-support solutions to matrix models in the planar limit ^Oj- A natural question 
which immediately arises in this context is whether these functions - the quasiclassical tau-functions, 
determined by multi-support solutions to matrix models, satisfy the WDVV equations? In the case of 
positive answer this is rather important, since multi-support solutions to the matrix models can play 
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a role of "bridge" between topological string theories and Seiberg-Witten theories which give rise 
to two different classes of solutions to the WDVV equations (see, e.g., ^21 and [T^ ITl] ). 

This question was already addressed in 15 , where it was shown that the multicut solution to 
one-matrix model satisfies the WDVV equations. However, this was verified only perturbatively and, 
what is even more important, for a particular non- canonical! (and rather strange) choice of variables. 

In this paper, we demonstrate that the quasiclassical tau-function of the multi-support solution 
satisfies the WDVV equations as a function of canonical variables identified with the periods and 
residues of the generating meromorphic one- form dS |16j . An exact proof of this statement consists 
of two steps. The first, most difficult step is to find the residue formula for the third derivatives ^ 
of the matrix model free energy. Then, using an associative algebra, we immediately prove that free 
energy of multi-support solution satisfies WDVV equations, upon the number of independent variables 
is fixed to be equal to the number of critical points in the residue formula. 

In sect. 121 we define the free energy of the multi-support matrix model in terms of the quasiclassical 
tau-function L16_ along the line of jl7l 1181 I19j . In sect. |51 we derive the residue formula for the 
third derivatives of the quasiclassical tau-function for the variables associated both with the periods 
{S} = {Si} and residues {t} = {ti} of the generating differential dS. In sect. 01 we prove that the 
free energy of the multisupport solution J^{T) solves the WDVV equations as a function of the 
full set of variables {T} = {S,t} whose total number should be fixed to be equal to the number of 
critical points in the residue formula for the third derivatives 0. In sect. El we verify this statement 
explicitly for the first nontrivial case where the total number of variables is equal to four. ^ Finally 
we present several concluding remarks and discuss possible generalizations. 

We restrict ourselves by the M = 1 supersymmetric theories without fiavours originally considered 
in 1^. The results can be easily generalized. Note that the literature on the subject is already quite 
extensive [201 E!, and different interesting developments of the issues discussed in this paper can be 
immediately obtained. 

2 Tau-function of multi-support matrix model 

We first exactly define what we call below the multi-support free energy of the matrix model in the 
planar limit. We are mostly doing with one- matrix integrals ^ of the form 

Z = l dcI>e^TrW^(*) (3) 

where the potential W^(<^) is supposed to be a polynomial of a degree (n -|- 1). The free energy in the 
planar limit of ^ can be defined as the first term in the expansion 

oo 

Fit, to) = lim(logZ(t)) = ^?i23-2^^(^^^^) 

9=0 

implying — > cxd, ^ ^ with Nh = to being fixed. In what follows we are only interested in the first 
term of this expansion, i<o(t,to). In fact, we deal with another quantity, J-'{t,to,S), where Si = hNi, 
J2 Si = to are extra variables - the filling numbers of (metastable) vacua. In order to get from this 
quantity FQ{t,tQ), one needs to minimize the free energy with respect to the filling numbers, = 0. 
However, one can still preserve Si as free parameters introducing the "chemical potentials". 

The origin of the new variables Si becomes rather transparent after one says that instead of direct 
computation of Q this problem is replaced by the saddle point approximation - finding the extremum 
of the functional Fq [p{X)] oc / Wp — J J /o(Ai) log |Ai — A2I /o(A2) -|- Hq (/ p — to) where Hq is just a 

^Let us point out that the WDVV equations ^ are nontrivial only for the functions of at least three independent 
variables. However, as we see below, the structure of residue formula for the matrix model free energy requires the 
minimal number of independent variables to be at least four! From this point of view, the origin of the "experimental 
observation" of 15 valid for a function of three variables still remains unclear to us. 

^The generalization to the two-matrix case |19| is rather straightforward and will be discussed in the last section. 
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Figure 1: Cuts in the A- or "eigenvalue" plane for the planar limit of 1-matrix model. The eigenvalues 
are supposed to be located "on" the cuts. The distribution of eigenvalues is governed by the period 
integrals S = p{X)dX along the corresponding cycles and the dependence of partition function on 
"distributions" Si is given by the quasiclassical tau-function ^ = ^ p{X)dX. 



Lagrange multiplier to fix the total normalization of the eigenvalue density. This latter condition 
means the saddle point equation is non-trivial only on the support of p. For one matrix model, this 
support can be presented as a set {Di} of cuts in complex eigenvalue plane, see fig. ^ Then, one 
should add to this functional the term (^Id P ~ '^«) ' '^^lich via Lagrange multipliers, controls the 

filling numbers at each cut, i.e. to consider 

(t, to, S) oc Jwp- J J p{Xi) log |Ai - A2I p{X2) +Uo(^J p- to) + E [J^P- (5) 

Then, the extra variables appear due to the extra information hidden in © compare to © - the 
structure of nontrivial eigenvalue supports. It is well-known that at "critical" densities §^ = 0, © is 
a (logarithm) of quasiclassical tau-function (see, e.g., [THIITU] !. 

In principle, in order to compare J- with the matrix model quantity Fq, one needs to put further 
restrictions to get rid of the metastable vacua. This would lead to shrinking part of the cuts into 
the double points (see discussions of these issues, say, in fTHI). Here we would forget this issue and 
consider smooth curves © with only two marked points at infinities on two A-sheets of the curve (0 . 

Note also that in ^ one can make two different natural choices for the set of new independent 
variables: the first choice corresponds to independent filling of all cuts, then to = "^Si, while the 
second choice corresponds to choosing as independent to and all Si except of corresponding to one of 
the cuts (that is why the corresponding sum in © is denoted as Y^'). These choices are related by the 
linear change of variables which does not influence the WDVV equations (see EOl ) • The first choice 
is more "symmetric" while the second one corresponds to the canonical choice of variables in the sense 
of jl6j or to the homology basis on smooth curve © with added marked points at two infinities. We 
use both of them below depending on convenience. 

The complex curve of one-matrix model "comes from" the loop equations (see, for example, |22j ) 
and can be written in the form 

y' = W'iXf + fiX)^RiX) (6) 
with the matrix model potential (jS)) parameterized as 

T^(A) = ^t;+iV+i (7) 

l>0 

or 

n 

W'iX)=Y,il + l)ti+,X' (8) 
1=0 
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being the polynomial of n-th degree in our conventions. The coefficients of the function 

n-l 

/(A) = Y: fkX' (9) 

fc=0 

are related to the extra data (the filling numbers) of the multicut solution. The eigenvalue density 
p{X) is the imaginary part of y(A) and vanishes outside the cuts. Therefore, the eigenvalue distribution 
S can be fixed by the periods 



S, = i dS (10) 

of the generating differential 

dS = ydX (11) 

taken around the eigenvalue supports to be identified (except for one of the supports) with the canonical 
A-cycles. Then 

ddS 

"^■^^ (12) 
dujj = 6ij 

when the derivatives are taken at fixed coefficients {ti} of the potential ©. One can show that the 
Lagrangian multipliers in ^ are given by integrals of the same generating differential (jllj) over the 
dual contours (see fig. ^ 

Ui= (f dS (13) 

Jb, 

To the set of parameters (|1()|) one should also add ^ the total number of eigenvalues Nh = to 

resoo (dS) = {"-^ ^ to (14) 
2(n+l)t„+i 

and the parameters of the potential 0, (jHI) which can be equivalently written as 



tk = -^resoo (a ^dS 
k = 1, . . . ,n 



(15) 



while the leading term coefficient t^+i is supposed to be fixed (we will discuss this issue in detail 
below). Then 

d^o = = [n + l)tn+i + oZ^^T (16) 

and the dependence of {fk} with /c = 0, 1, . . . , n — 2 on to is fixed by the condition 

n + 1 t„+i + ^ ^ = 17 

which for i = 1, . . . , n — 1 gives exactly n — 1 relations on the derivatives of /o, /i, . . . , fn-2 w.r.t. to. 
The bipole differential (jlGj) can be also presented as 

^By the oo-point in what follows we call for short the point co+ or A = oo on the "upper" sheet o f hyperelliptic 
Riemann surface Q corresponding to the positive sign of the square root, i.e. to y — +a/W^'(A)^ + /(A). 
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where E{P, P') is the Prime form. Differential ()18|1 obviously obey the properties 



veSood^Q = — resoo_d^^o = 1 

(19) 



<f dr^o = 0, i = l,...,n — 1 

For the derivatives w.r.t. parameters of the potential (|15)) . one gets 

ddS _ W'{X)k\''-^d\ 1 dfj \^d\ 



obeying 



dtk y 2^^dtk y 

A, J A, y dtk Ja, y 

and this is again a system of linear equations on To complete the setup one should also add to 
(|13() the following formulas ^: 

no= / dS (22) 



OO- 

n-1 



(we again remind that, instead of t^, the parameter Sn = tQ — Yl^=i Si can be used equivalently) and 

Vk = resoo (x^dS) , k > {) (23) 

On genus g = n—1 smooth Riemann surface ©, there are 2g = 2n—2 independent noncontractable 
contours which can be split into the so-called A = {^i} and B = {-Bj}, i = 1, . . . ,g, cycles with the 
intersection form Ai o Bj = 6ij. The canonical holomorphic differentials p2j) are normalized to the 
A-cycles, and their integrals along the B-cycles give the period matrix, 



duJi = Tij (24) 

To check integrability of ()13() and (|22() one needs to verify the symmetricity of the second derivatives. 
For the part related with the derivatives only w.r.t. the variables (jlUj) . this is just a symmetricity of 
the period matrix of © and follows from the Riemann bilinear relations for the canonical holomorphic 
differentials ((T^ 



= / duJi A dujj = ^ I i duJi j> dujj — j> dujj j> duJi J = Tij — Tji (25) 



Analogously 



= / duJi A dQo = ( (p doji f dViQ — f dVi^ f duJi ) + 
Jt. ^ \JAk JBk JAk JBk 

+reSoo(d^2) / — resoo(df^o) / duJi = d£l^ — \ duJi 

Joo— J CO— ^ Bi J CO — 



(26) 



Formula H25|) means that 



Qg- - T,^ - T,, - — (27) 



* Naively understood the integral in 12211 is divergent and should be supplemented by some proper regularization. 
In what follows we ignore this subtlety since it does not influence the residue formulas for the third derivatives, those 
one really needs for the WDVV equations The simplest way to avoid these complications is to think of the pair of 
marked points oo and oo_ as of degenerate handle; then the residue (1141 comes from degeneration of the extra A-period, 
while the integral 1221 from degeneration of the extra -B-period. 
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Figure 2: Cut Riemann surface with boundary dT,. The integral over the boundary can be divided 
into several pieces (see formula (|32])). In the process of computation we use the fact that the boundary 
values of Abelian integrals on two copies of the cut differ by period integral of the corresponding 
differential dujj over the dual cycle. 



while from (|26j) one gets 

dUj _dUo 

~dt^-dS~ ^^^^ 
This allows one to introduce the function J^{T) = J-{S,tQ,t) such that 

wr""^^ ^r""'^ drr'' ^''^ 

The integrability of the last relation can be checked similar to l\'27\i. H28|) with the help of Riemann 
bilinear relations involving the Abelian integrals Qk = dO.^, for example (cf. e.g. with [S^, where 
similar relations were used for the quasiclassical tau- function of the Seiberg-Witten theory): 



resn 



{0,kduJi) = j> ^ flkdcoi = ^ (J^ ^'IduJi - J^Qf^ duJi) - ^ (J^ ^^duJi - duoi 

= 2_, ( f dQ,k f duJi — f dUk f duJi I = f dO.^ 
, \JBi JAi JAi Jbi J Jb, 



(30) 

where dT is the cut Riemann surface © (see fig. \^ , and in the last equality we used (|21|1 . 



3 Residue formula 

3.1 Holomorphic differentials 

Let us now derive the formulas for the third derivatives of following the way proposed by I.Krichever 
jl6l I25j . We first note that the derivatives of the elements of period matrix (in this section, for 
simplicity, we set the coefficients of potential (jSI) to be fixed) can be expressed through the integral 
over the "boundary" dT of cut Riemann surface T (see fig. |2) 

= dkTij = I dkdLOi = I ujjdkdLOi (31) 
dSk Jb, JdT. 



6 



where ujj = dujj are the Abehan integrals, whose values on two copies of cycles on the cut Riemann 
surface (see also fig. are denoted below as uj^ . Indeed, the computation of the r.h.s. of (|^T|) gives 

/ ujjdkdui = X! ( / ^fdkduji - / ujjdkduji ) - X! ( / ^f^kduji - / uJdkdiVi \ = 

J dYli ^ V-^ Bi J Bi / ^ \^ A.I J A.I J 

= ? i a - ? L a = (32) 

= y^A f dujj] (f) dkdui - V ( <P dwj] f dkduJi = dkTij 
^ \JAi J JBi I \JBi J JAi dwj=5,j 

One can now rewrite as 

dkTij = j^^ LVjdkdwi = - j^^ dkujjduji = ^ resdA=o {dkUJjdoJi) (33) 

where the sum is taken over all residues of the integrand, i.e. over all residues of dkUJj since dwj are 
holomorphic. In order to investigate these singularities and understand the last equality in (|33|) . we 
discuss first how to take derivatives dk w.r.t. moduli, or introduce the corresponding connection. 

To this end, we introduce a covariantly constant function - the hyperelliptic co-ordinate A, i.e. 
such a connection that dkX = 0. Roughly speaking, the role of covariantly constant function can be 
played by one of two co-ordinates - in the simplest possible description of complex curve by a single 
equation on two complex variables. Then, using this equation, one may express the other co-ordinate 
as a function of A and moduli. Any Abelian integral loj can be then expressed in terms of A, and 
in the vicinity of critical points {Aq} where dX = (for a general (non-singular) curve this is always 
true) we get an expansion 

UJj{X) = UJja + CjaV^- K + ■■■ (34) 

A — >^Aq- 

whose derivative 

dkujj = dkUJj I ^^^^^^^ = - ^^^'^ dk \a + regular (35) 

gives first order poles at A = Aq, up to regular terms which do not contribute to (jSSI). The exact 
coefficient in can be computed for A related with the generating differential dS = ydX. Then, 
using 

y(A) = r^VA- A„ + ... (36) 



where r„ = JUf^jLai^a - A/j) or 



together with 



and 



dy = — " dX + regular (38) 

2^/ X — Aq 



rf^. = ^^7==^rfA + ... (39) 



and, following from and H12|) 



diok = dkdS = " ^ " dX + regular (40) 

2 V A — Aq, 



one finally gets for (|33|) 
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In hyperelliptic situation, the derivation presented above is equivalent to using the Fay formula 



dX. 



BT- ■ 

-Q:i{K)(^^{K) (42) 



where Q^iiXa) = ^\^\ is the "value" of canonical differential at a critical point. 



3.2 Meromorphic differentials 

Almost in the same way the residue formula can be derived for the meromorphic differentials pflj) . 
One gets 

1^ = resoo {\^dS\ , > (43) 

therefore 



a. ^^'Soo = reSoo ((nfe)+dn„) (44) 

where (J^a;)^ is the singular part of the integral of 1-form d£L\j,. Further 

— — reSoo {^(KXri\ = reSoo f^A^— - — = — reSoo { {d^k) + -?r—- \ = — reSoo ( dQk-?r—\ (45) 
The last expression can be rewritten as 

— reSoo ( dQh ^ " ^ = / ( dClk^r—-^ +y^i"es\ ( dClu ^^ =y^i'esA ( dVtu ^^ (46) 
since (^^fc^) = due to dO„ = 0, (cf. with 

= i: i ii "".) - E £, (i ^".) il"". = 



Now, as in the holomorphic case one takes 



and, therefore 

d _ d 



X=const 2\/ A — Aq, 9t 

Then, using (1211), (jSZI) and (jHEl) together with 



and the relation, following from (|2fl|) . (|.S7j) 



(47) 



On(A) = r^na + 7n«VA - Aq + . . . (48) 

A — >Aq, 



+ regular (49) 



^% = ^^/g^ + (50) 
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one gets for 



resA, dJ],. — - — = resA, 



(iA " V dXdy 

In a similar way, one proves the residue formula for the mixed derivatives, so that we finally conclude 
d^T ^ fdHidHjdHK\ V- (^i(t)jcpK 



dTidTjdTx 



^ f dHjdHjdHK\ ^ f 4>i4>j(l>K r \ 



. . . (53) 

for the whole set {Tj} = {tk,tQ, Si} and {dHj} = {drifc, dOo, duJi}. In formula (|53jl we have introduced 
the meromorphic functions 

for any (meromorphic or holomorphic) differential dHj = Hi{X)^ and -R(A) = W^'(A)^ + /(A). The 
derivation of the residue formula for the set of parameters including to corresponding to the third-kind 
Abelian differential H16|) can be performed in a similar way. 

4 Proof of WDVV 

Having the residue formula , the proof of the WDVV equations is reduced to solving the system 
of linear equations [7| I26[ IT^. which requires only fulfilling the two conditions: 

• The "matching" condition 

#(/) = #(«) (55) 

and 

• nondegeneracy of the matrix built from ()54l) : 

det||(/>/(A„)|| ^0 (56) 

la 



Under these conditions, the structure constants C^j of the associative algebra 

)f iCl)M 

(57) 



(Ci)'} ^ Cfj 



responsible for the WDVV equations can be found from the system of linear equations 

(t)i{K)(t>j{K) = E CfjcpxiXa), V Xa (58) 



K 

with the solution 



Cfj = E '/'/(Aa)0j(Aa) {<t>K{Xa)r' (59) 

a 

To make it as general, as in Pl IH 01 , one may consider an associative isomorphic algebra (again V Aq,) 

4>l{Xa)(t>j{K) = E CfjiO^KiXa) ■ e(Aa) (60) 



K 
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which instead of (fS^ leads to 

Cf.(0 = E^%W^(^-(A.)r (61) 
The rest of the proof is consistency of this algebra with relation 

:FijK = Y.Cu{i)mL{i) (62) 

L 

with 

riKLii) = iA^KLA (63) 

A 

expressing structure constants in terms of the third derivatives and, thus, leading to (^. It is easy to 
see that (|62|) is satisfied if Tkla are given by residue formula (jSHJ. 
Indeed, requiring only matching #(a) = #(/), one gets 

T.cfAi)mL{K) = E '^'^^"af^"^ • ('^^(^-))"' • ^K{x,)M\,m,n (64) 

K K,a,f3 ^'^^"^ 

and finally 

E CfAOVKLiO = E '^'^^"a^/^"V L(Aa)e(Aa)ra = J2^aMK)MK)^LiK) = TljL (65) 



K 



Hence, for the proof of one has to adjust the number of parameters {T/} according to (|55|) . The 
number of critical points #(a) = 2n since dA = for = R{X) = 0. Thus, one have to take a 
codimension one subspace in the space of all parameters {T/}, a natural choice will be to fix the eldest 
coefficient of 0. Then the total number of parameters #(/), including the periods S, residue to and 
the rest of the coefficients of the potential will he g + 1 + n = (n — l) + l + n = 2n, i.e. exactly equal 
to #{a). In sect. El we present an explicit check of the WDVV equations for this choice, using the 
expansion of free energy computed in ^ |^ . 

Note that equations (|6fl|) are basically equivalent to the algebra of forms (or differentials) considered 
in [HI 15 . In this particular case one may take the basis of 1-differentials duJi , dQ^ with multiplication 
given by usual (not wedge!) multiplication modulo dS = ydX. Then, one can either directly check 
that the algebra with this multiplication is associative (similar to how it was done in I^EI), or, using 
hyperelliptic parameterization, remove the factor ^ in order to reduce the algebra to the ring of 
polynomials with multiplication modulo the polynomial ideal = H^'^(A) + /(A), which is obviously 
associative. 

In the proof of the WDVV equations we used in this section, the associativity (|57|) is absolutely 
evident, being associativity of the usual multiplication, and the main point to check was to derive 
(|62j) - (|63|) . When using instead the algebra of differentials, the main non-trivial point is to check its 
associativity, while the relations (|62() - H63|) appear as even more transparent than above corollary of 
the residue formula. 



5 Explicit check of the WDVV equations 

To convince ourselves that the general proof indeed works and to get some further insights, in this 
section we consider the explicit check of the WDVV equations (0) perturbatively. To do this, let us 
take the perturbative expansion of the prepotential ()29|) at small /Sj's (here we take the symmetric 
choice of variables with i = 1, . . . ,n = g + 1), keeping the coefficients of the matrix model potential t's 
arbitrary, and see if this perturbative expansion of satisfies the WDVV equations order by order. 
A general procedure of getting such perturbative expansion was proposed in j^. 
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We are going to check here only the simplest non-trivial case of cubic matrix model potential, 
which is reformulated in terms of elliptic curve © with four branch points. According to (jSSj) the 
corresponding solution to WDVV equations should depend exactly on four independent variables, 
and we are choosing them to consist of two filling numbers {Si and ^2) and two coefficients of the 
potential (ti and 12). The perturbative expansion for this case was constructed in (9j up to the fifth 
order in S'j's and was later discussed in many other places (see, e.g., It reads 

= -SiWi^ii) - S2W{fi2) + ^Sf log (^^) + ^Sl log (^^) + 2S1S2 log(A)+ 

+^ [Isf - + 5S,Sl - 5S!S2) + ^ (^Isf + - fSiSl - |5?52 + 595?5i) + (66) 

,1 /^56 5 56 5 871 4 871 4 . 2636 3 2 ^^Sf) 2 c.3\ , qSn 
A9 VT ~ "3" 3~ ^ -^•^2'^i + — ^'^i'^2 3~ 1 2 ) + ) 

Here the matrix model potential ((7)) is fixed to be W{X) = A'^/3 + t2X'^ /2 + tiX; and fii, ^2 are the 
roots of the equation W'{X) = 0, i.e. 

^2 = -\t2-\^tl-^t, (67) 
A = ^1 - /i2 = \lt2- 4*1 

The perturbative expansion of free energy (|66|) is symmetric with respect to simultaneous transforma- 
tion S\ <-> —52 and t2 <-> — ^2- However, since the whole prepotential but its linear part depends only 
on the transformation of ^2 is only essential for the linear term. 

Note that in the context of supersymmetric field theories such formulas usually depend on an 
additional parameter ^qcd associated with the field theory scale However, this parameter here 
emerges only as regularization (see footnote^ and generally can be omitted within the matrix model 
approach. Moreover, since the terms depending on Kqcd are at most quadratic in Si and ti (they are 
logarithmic or polynomial in kqcD^ with the latter proportional to PF(Aqc'£))), they do not contribute 
to the WDVV equations and just renormalize the second derivatives of the prepotential (or period 
matrix of the curve similar to © playing the role of the set of effective couplings in the field theory 
at low energies (cf. with _11J). 

There is one subtlety with the perturbative prepotential - when calculating it the authors 
of were interested only in the terms depending on Si. There could be, in principle, some terms 
dependent only on tj's. This does not, however, happen in our case. Indeed, these terms would survive 
in the limit of all Si = 0, i.e. when /(A) = in (jB)) and y = W'{\) is just a polynomial. Then from 
(|23j) it follows that Vk = = in this limit, hence we miss at most linear in ti terms. 

In contrast to the weak coupling expansion of the Seiberg-Witten theory when logarithmic terms 
do satisfy the WDVV equations themselves exactly (see 0)^, some technical complication with pertur- 
bative check of whether the prepotential satisfies the WDVV equations emerges due to different 
orders of magnitude of different matrix elements of the matrix of third derivatives of the prepotential 
((21). Say, matrix {^12) jk has all matrix elements except for (^ta)*^ t^. the order 0(1), while the 
latter ones of the order of 0{Si). (We definitely put all Sj's of the same order of magnitude.) 

Therefore, a careful calculation in every matrix element is required to determine the order of 
the WDVV equations where the perturbative expansion contributes to. In practice, one suffices to 
rescale all 5i's with a scale parameter A. Then, all non-trivial matrix elements of same matrix WDVV 
equation (i.e. for concrete UK) are of the same order in A (choosing other indices UK in changes 

^Note that expansion 15511 that goes over positive powers of variables, Si's is rather similar not to the weak coupling 
expansion of the Seiberg-Witten theory but to its strong coupling expansion (see, for example, |29|). In the latter case, 
the WDVV equations trivially hold due to the duality argument |30|. 
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this order). Say, for the choice {I,J,K) = {t2,Si,ti) the leading order of the equation(^ is and 
trivially vanishes. It is contributed only by the linear and logarithmic (in S'j's) terms of ()66() . The 
cubic terms of the prepotential contribute to A^-order, the quartic terms to A^-order etc. Using the 
program MAPLE, we have checked that the WDVV equations are satisfied with the prepotential 
l|66|) up to the fifth order in 5j. 

Now, one can consider further non-trivial examples with matrix model potentials being polynomials 
of higher degree. One technical problem now is that W{^i) is getting more and more involved function 
of tj's, since //j's become the roots of polynomial equations of increasing order. Therefore, we lose 
just a principal possibility of analytic analyzing W{ni) if VF(A) is a polynomial of order 6 and higher. 
Nevertheless, one can still look at different limiting regimes with some of roots of W'{X) approaching 
close to each other. 

6 Conclusion 

In this paper we have proven that the (generalized) free energy of the matrix model in planar limit 
satisfies the WDVV equations . The proof is based on the residue formula for the quasiclassical tau- 
functions introduced in ^S] which was derived above for the whole basis of the first-kind or holomorphic 
(|12jl . second-kind (|20|) and third-kind (|16j) meromorphic Abelian differentials. The proof of the WDVV 
equations based on the residue formula requires the matching condition (or associativity of the 
algebra of differentials) , and the first nontrivial example of our statement for the case of four variables 
we have checked explicitly for its perturbative expansion H66() . 

The matching condition (|55|) requires for the quasiclassical tau-function defined by © and (|llj) 
to satisfy the WDVV equations (pQ) as a function of the set of parameters involving necessarily the 
coefficients of the matrix model potential and corresponding to the meromorphic differentials. In other 
words, we cannot restrict ourselves for the set of only holomorphic differentials, like it happens for the 
prepotentials of the Seiberg-Witten models [2], associated as well in the sense of sect.|21to integrable 
systems |27| (see also [2E1 and references therein). This is similar to the case of softly broken = 4 
Seiberg-Witten theory, where the WDVV equations are not satisfied by the prepotential, being a 
function of the Seiberg-Witten periods corresponding to the holomorphic differentials only In 
that case one should also necessarily add meromorphic differentials, but any physical sense of the 
corresponding parameters, in contrast to the situation considered in this paper, remains yet unclear. 

Our statement can be generalized straightforwardly to the case of the two-matrix model and the 
corresponding non-hyperelliptic curve JHI- It was already suggested in [TP that the WDVV equations 
should hold for the quasiclassical tau-function of the two-matrix model, whose construction is based 
on the non-hyperelliptic curve, and similar to (lllj) meromorphic differential. One can use our line 
of reasoning of sect. |31 and sect. 21 as a general proof in the two-matrix case either, since our proof 
does not require any specific of the curve © like existence of the hyperelliptic parameterization. The 
perturbative expansion of the multisupport solution of the two-matrix model is not known yet (it can 
be computed, say, using the diagrammatic expansion proposed in |19j ) and it would be interesting 
to calculate the free energy of the two matrix model explicitly and repeat for this case the explicit 
calculation of sect. [51 

There is another, very important direction where one can easily extend the consideration of the 
present paper. In there was shown how one can introduce more Whitham times, preserving the 
same number of moduli Si. It is based on using potential Ty(A) of higher degree and fixing some of the 
coefficients of /(A) by the double point conditions. This procedure allows one to make the number of 
Whitham times and moduli independent in the WDVV equations. It would be interesting to check if 
the WDVV equations are satisfied in this case, and to construct the perturbative expansion for such 
a case. 

Finally, let us point out that the fact that free energy of multi-support solution of matrix model 
does satisfy the WDVV equations can be thought of as a first step to write down the generalization 
of the (quasiclassical) Hirota equations which are known up to now only for the one-support solutions. 
A link between the WDVV equations and (generalized) Hirota relations [3 ESI together with the 
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statement of our paper allows one to believe that the counterpart of the Hirota relations for nontrivial 
quasiclassical tau-functions can be written in an explicit form. 
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